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Abstract 

We study the chiral behavior of the nucleon and A-isobar masses within a manifestly covariant 
chiral effective-field theory, consistent with the analyticity principle. We compute the ttN and ttA 
one-loop contributions to the mass and field-normalization constant, and find that they can be de- 
scribed in terms of universal relativistic loop functions, multiplied by appropriate spin, isospin and 
coupling constants. We show that these relativistic one-loop corrections, when properly renormal- 
ized, obey the chiral power-counting and vanish in the chiral limit. The results including only the 
7riV-loop corrections compare favorably with the lattice QCD data for the pion-mass dependence 
of the nucleon and A masses, while inclusion of the ttA loops tends to spoil this agreement. 
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The nucleoli mass (Mat ~ 940 MeV) is much larger than the sum of the masses of its 
constituents (3m q ~ 20 MeV), hence almost all of it is generated by the strong interaction 
among the quarks. An exact description of this phenomenon has not yet been derived from 
QCD, however, tremendous progress has been achieved in the numerical computation of the 
nucleon mass in lattice QCD 0,0]. One of the main limitations of the lattice QCD studies is 
that the finite lattice size restricts the value of quark masses from below, and thus the quarks 
in present lattice studies are much heavier than in reality. It became a common practice 
to perform lattice calculations for different values of quark masses and then extrapolate the 
results to the physical point. 

The extrapolation in the quark mass is not straightforward, because the non-analytic 
dependencies, such as ^/m q and \nm q , are shown to be important as one approaches the 
small physical value of m q . Therefore naive extrapolations often fail, while spectacular 
non-analytic effects are found in a number of different quantities, see e.g., Refs. H 0, E| . 
Fortunately, it is known how to compute these non-analytic terms in chiral effective field 
theory (ChEFT) — a low-energy effective field theory of QCD. For recent examples of such 
calculations for the nucleon and other baryon masses see, e.g., Refs. @, 0, i, B H E2 



In this Letter we present a new calculation of the quark-mass dependence of the nucleon 
and A- isobar masses in the framework of relativistic ChEFT, with the emphasis on the 
analyticity constraint. 

In the ChEFT the interaction is mediated by pions, which are the Goldstone bosons of 
the spontaneously broken chiral symmetry of QCD. The explicit-chiral-symmetry breaking 
terms, represented by the pion and quark masses in ChEFT and QCD, respectively, are 
related via the Gell-Mann-Oaks-Renner relation: f%m^ = — (qq)2m q , where (qq) ~ —(230 
MeV) 3 represents the value of the quark condensate. Lattice calculations confirm this re- 
lation for a very broad range of quark masses |l2|. Thus, the quark- mass dependence of 
quantities in QCD can be translated to the pion-mass dependence of these quantities in 
ChEFT and vice versa. 

As the strength of the Goldstone boson interactions is proportional to their energy, at 
sufficiently low energies a convergent perturbative expansion in ChEFT is possible. However, 
most of the lattice results are presently obtained for pion masses above 300 MeV where the 
chiral expansion is not expected to converge well. Therefore one resorts to methods where the 
leading non-analytic ChEFT results are combined with more phenomenological techniques 
such that the resulting approach has a wider range of applicability 0, albeit lesser predictive 



power. 

Recently it has often been arg ued EI ELI [03 that the manifestly relativistic ChEFT cal- 



culations have, in some cases, better convergence than their heavy-baryon (semi-relativistic) 
counterparts. This implies that the convergence of the ChEFT expansion is improved by 
a resummation of nominally higher-order terms which are relativistic corrections to the 
leading non-analytic terms. One can thus improve on the convergence of the chiral expan- 
sion without loss of predictive power, or in plain words, without introducing additional free 
parameters. 

The original formulation of chiral perturbation theory with nucleons had been relativis- 
tic Esl> but was claimed to violate the chiral power counting. The so-called heavy-baryon 
chiral perturbation theory, which treats nucleons semi-relativistically, was developed to cure 
the power-counting problem E3|' an d a l°t °f work has been done since in this direction. 
More recently, Becher and Leutwyler l| proposed a manifestly Lorentz-invariant formu- 
lation supplemented with so-called infrared regularization (IR) of loops in which the chiral 
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power-counting is manifest. At about the same time it was realized |18( that power-counting 
can be maintained in a relativistic formalism without the IR or the heavy-baryon expan- 
sions. The original, straightforward formulation complies with chiral power-counting if 
appropriate renormalizations of available low-energy constants are done. 

Power-counting issues apart, the original relativistic formulation has a particular advan- 
tage over the IR scheme in that it preserves analyticity of the loop contributions . The 
IR procedure spoils the analyticity by introducing unphysical cuts in the complex energy 
plane. In our work we therefore prefer to use a manifestly Lorentz-covariant formulation of 
ChEFT, supplemented with appropriate renormalizations ^H, rather than infrared regular- 
izations 13], to maintain power counting. 

We begin with defining the effective chiral Lagrangian. Writing here only the first-order 
terms involving the isovector pseudoscalar pion field 7r a , the spin- 1/2 isospin-1/2 nucleon 
field N and spin-3/2 isospin-3/2 field ^ of the A-isobar we have (in the conventions of 
Appendix A): 

= N( % @-M N )N + -^-N 1 ^ l5 T a (d ll N)d u n a 

+ ^^ a id a - M Al n $ v + e^ aX ^J a {d a ^) d X 7r a (1) 

where Ma — 1232 MeV is the A-isobar mass, f n ~ 92.4 MeV is the pion decay constant, 
qa — 1.267 is the axial coupling of the nucleon, while Ka and Ha represent the lowest 
order 7riVA and 7rAA couplings, respectively. In the large-Ac limit they are related to 
qa as Yla = (3/v / 2)fi' J 4, Ha = 9/5<^. The isospin factors enter through the Pauli matrices 
r a , the isospin- l/2-to-3/2 transition matrices T a , and the isospin-3/2 matrices T a , with 
normalizations, r a r a = 3, T a ^T a = 2, T a T a = 5/3, where summation over a (= 1,2,3) is 
understood. 

To study the chiral behavior of the nucleon and A masses we introduce a counter-term 
Lagrangian containing the corresponding quantities in the chiral limit (m n — > 0): 

45. = (M N -M^)NN+[zi%-l]N(i^-M^)N 

+ (M A - Mf) %-T1> v + [Zfl - 1] fyW^id,, - Mi 0) 7 ^) ^ , 

Cfl = 4c 1N m 2 n NN - 4d 1N mlN{i$ - M^)N (2) 
+ 4 c 1A ml - 4 d 1A mffi fi {i a ' a id a - Nlf 1 ^) ^, 

where and represent the chiral-limit value of the masses and the field- 

renormalization constants, respectively. 

Our choice of the chiral Lagrangian is different from the ones previously used in the 
literature (e.g., [l(], 11, Hf) in two important aspects: 



(i) The ttNN coupling differs from the usual pseudovector coupling: Jj- N (<^7r a )7 5 r a N , 
which is standardly used at this order. The difference between this and our nNN 
coupling is of higher order as can easily be shown by using partial integration and the 
Dirac equation for the nucleon field. Nonetheless, our choice simplifies the calculation 
and, most importantly, allows us to write down the results for the nucleon and the A 
in the same form, see Eq. (|T7|) below. 
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(ii) The couplings of the spin-3/2 field are invariant under a gauge transformation: 



-> Vv + fye, (3) 

with e a spinor field. This requirement is called for by the consistency with the free 
spin-3/2 field theory (l9| . which is formulated such that the number of spin degrees of 
freedom is constrained to the physical number, see Refs. [13, Hl[ for details. 



Both of these points are crucial for the consistency and elegance of this calculation. 

Point (ii), in particular, allows us to use simpler forms for the spin-3/2 propagator. 
Indeed, as can be read off Eq. ((TJ), the propagator of the massive spin-3/2 field is the inverse 
of the free-field operator: 

(S~ l )a(3(p) = la/3^P^ -m^ap, (4) 

where p = id, and m denotes the mass. However, using the gauge symmetry under Q and 
hence the spin-3/2 constraints: d • ip = = 7 • ip, one can obtain other, equivalent, forms of 
the propagator |22|. One can, for example, derive the following gauge-fixing term: 

£ g .f. = -*C (P ■ $ 7 • V> - $ • 7 9 ■ if>) , (5) 

with the gauge-fixing parameter (, a real number. Upon adding this term, the free-field 
operator Eq. (jlj) becomes: 

(S^apip) = (p'-m) 7 Q(8 + (1 + 0(laP/3 ~ 1/SPa) 

= la/3 {j/)-m)-{l- 0(laPf3 ~ IfBPa) , (6) 



and it is not difficult to find its inverse: 

5 Q/3 (p) 



i2 „2 



771 — p 

Some simple gauges are: 



3(C p — m z ) i{Q P — wn ) 



(7) 



C=l: (8) 
C = -l: ^)=(^-hY)^, (9) 

C = oo : s^(p) = ^±^- V tm^ {p)) (10 ) 

m z — p 2 

V^{p) = | (> - + J^p, (11) 

is the covariant spin-3/2 projection operator. Obviously, ( = corresponds with the usual 
Rarita-Schwinger propagator. It is interesting to observe that for ( ^ the propagator has 
a smooth massless limit. We would like to stress that our results are independent of the 
gauge-fixing parameter, because all the spin-3/2 couplings used here are symmetric with 
respect to the gauge transformation (jBJ). 



where 
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FIG. 1: The nucleoli and A self-energy contributions considered in this work. Double lines represent 
the A propagators. 
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TABLE I: The coefficient Cbb 1 entering the ttN- and 7rA-loop contributions in the baryon mass 
formula Eq. Q17|). The rational numbers in the brackets represent the spin (= isospin) factors. 



In this spin-3/2 formalism the A self-energy takes a simple form: 



(12) 



where S{J>) has the spin-1/2 Lorentz structure. Thus, both nucleon and A-isobar self- 
energies can be expressed in the same Lorentz form, without complications of the lower-spin 
sector of the spin-3/2 theory considered in 24 . 

Furthermore, in explicit calculations we find that this form for the nucleon and the A can 
be written in a universal expression. Namely, the one-pion-loop contribution of a baryon B' 
to the self-energy of a baryon B, see Fig. ^ can generically be written as: 



c 



BB> 



B\ 



d A k 



3(2f 7T M B ) 2 i J (2vr) 4 k 2 - m\ (p - k) 2 



M 2 , 



[p 2 k 2 -(p-k) 2 ] 



(13) 



where Cbb' is given by the corresponding coupling constant squared, multiplied by the spin 
and isospin factors, see Table UJ To bring the spin-3/2 A-isobar contributions to this form, 
the identities listed in Appendix A are helpful. The similarity of the nucleon and A-isobar 
loop contributions, pointed out earlier in Ref. j2^|, is thus obtained here in the formalism of 
relativistic baryon ChEFT. 

To evaluate the loop integral we use the standard Feynman-parameter trick: (AB)^ 1 = 
f^dx [xA + (1 — x)B}~ 2 , and after the change of variable k — > k + (1 — x)p, obtain: 



C 



BB' 



3(2UMb) 



dx (xp 1 + Mb>)(p g 



2 iw 



jfp V )- 

i 



d A k 



kfi k v 



(2tt) 4 (k 2 - M 



2\2 



(14) 



mix + M§,(1 — x) — x(l 



with M 2 

dimensional regularization (for d — > 4r 



x)p 2 . The latter integral can be computed via 



d d k 



k/j,k u 



(2ir) d (k 2 - M 2 y 



2(4tt) 



4-d 



+ 7b - 1 - ln47r + \n(M 2 /A 2 



(15) 
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1 24^^» A = 


\ 0, A = 


M A 
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1 -96^^' A = 


I 2{^U)^ h A^ A > m ^ 

\ 0, A = 


3 25 rr2 

32vr/2 81 A 






TABLE II: The leading non-analytic contributions to the nucleon and A masses from the ttN- and 
7rA-loop. 



where the Euler constant 'Je = — r'(l) ~ 0.5772, and A is a renormalization scale. 

Writing the self-energy in general as E(j>) = <j(s) + {jf> — Mb)t(s), with s = p 2 , we find 
that 

*{*) = ~^JJ2^2 j/ x ( xM b + Mb,) M 2 [I. + HM 2 /s)] , (16a) 

r(s) = -^^^£dxxM a [l M + HM 2 /s)], (16b) 

with l s = —2/(4 — d) + 7_b — 1 — ln(47rA 2 /s). Obviously, a(M B ) contributes to the mass of 
baryon B, while r(M|) contributes to its field-renormalization constant (FRC), namely: 

M B = mW-^^-j^j^Cm-^Hj,^), (17a) 

Z 2B = 42 -id lB ml + -0L_^c BB .W(ft,*>5^), (17b) 

where functions V and given explicitly in Appendix B, represent the one- loop contri- 
butions with mj and m 2 terms subtracted. The latter terms are subtracted because they 
merely renormalize the available low-energy parameters, here Z% , ci, and d±. It is 

interesting to note that the loop function V is a relativistic analog (up to a constant factor) 
of the function W of Banerjee and Milana which represents the heavy-baryon results. 
For Mb = Mb 1 the loop functions simplify considerably: 

3 

V(fx, 0) = — {8 (1 - i/i 2 ) 5/2 arccos £ + £ [17 - 2^ 2 + (30 - 10/i 2 + /i 4 ) ln/i 2 ] ) ,(18a) 
3 1 2 8 j 

3 

Wfo 0) = ^ {2 (1 - i/i 2 ) 3/2 arccos £ + £ [13 - 2^ 2 + (18 - 8/i 2 + /1 4 ) ln/i 2 ] ) . (18b) 
3 1 2 8 J 

By studying the expansion of these functions near the chiral limit, see Appendix B, we 
find that the chiral expansion for the mass goes as: 

4 

M B = Mg ] - 4c 1B ml + XBBml + XBB'^lnm^ + 0{ml), (19) 

where A = Ma — Mjy > m v and the chiral coefficients can be read off Table |H] This 
expansion shows explicitly that the (renormalized) relativistic loop contributions vanish in 
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FIG. 2: (Color online) Pion-mass dependence of the nucleon mass. The dashed (blue) curve is 
the leading- nonanalytic m 3 result, whereas the solid (black) curve is the full relativistic 7riV-loop 
result, both for parameter values: = 0.883 GeV and c\n = —0.87 GeV -1 . The dotted (green) 

curve is the relativistic result for ttN + 7rA loops with M$ = 0.87 GeV and c\n = —1.1 GeV -1 . 
Upon adding to this result the m\ term, Eq. (|2^|) . with C2N = 3.0 GeV -3 , one obtains the dash- 
dotted (red) curve. The (red) squares are lattice results from the MILC Collaboration [2||. The 
star represents the physical mass value, which is used in the fits. 



the chiral limit, as they must. Also, since there are no contributions which, near the chiral 
limit, scale with positive powers of A, the introduction of counter-terms of such nature, as 
is done in is unnecessary, and would be excessive in this calculation. 

The expansion (J19)) also shows that the loop contributions obey the chiral power counting. 
For example, in the so-called 5-counting [26{, the ttN and ttA loop contributions to the nu- 
cleon mass count as p 3 and p 4 /A, respectively, which for small p ~ m n agrees with Eq. (|19j) . 
In the 5-counting the one-loop result Eq. (JT7J) thus represents a complete calculation to order 
p 4 /A ~ 5 7 . A full fourth (p 4 ) and p 5 /A order calculation of both ttN and 7rA loops requires 
calculation of diagrams in Fig. ^ with a ttBB' vertex from (two derivatives of the pion 
field) and the tadpole contributions. Such a calculation is a worthwhile topic for a future 
work. 

The nN contribution to the A self-energy has an imaginary part for m T < A, which gives 
rise to the A width. According to this calculation the width is given by [27]: 

fA = UMliiu)^ + MN? ~ ml?/2{A2 ~ mlf ' 2 - m 

The experimental value, i~A — 115 MeV, fixes Jia — 2.85, the value which we shall use in 
numerical calculations. Note also that this value is in a much better agreement with the 
large- Nc value Ha = 3gA/V% — 2.70, than with the SU(6)-relation value, Ha = 6gA\/2/5 — 
2.15. For the A axial coupling Ha, we use the SU(6) relation, which in this case coincides 
with the large-iV<7 relation: H A = (9/5)gA — 2.28. 

We are now in position to discuss the numerical results. Fig. E] displays the pion-mass 
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dependence of the nucleon mass, as given by Eq. (|17|) . The two low-energy constants M N 
and cin are related to reproduce the physical nucleon mass at the physical pion mass value. 
The only free parameter can then be adjusted to reproduce the lattice data, shown by 
the squares. Note that these lattice data are not corrected for finite volume effects, which 
are known to increase with decreasing m n , and have been estimated to reach 0.03 GeV 
for m 2 = 0.1 GeV 2 [9]. The solid curve in Fig. |2] shows the 7riV-loop contribution to the 
nucleon mass, with Mffl = 0.883 GeV and C\n = —0.87 GeV -1 . Thus, the relativistic 
calculation is able to describe the lattice results up to m 2 ~ 0.5 GeV 2 with only a single 
free parameter. For comparison, the dashed curve shows the corresponding leading non- 
analytic result [m 3 term in Eq. (jl9j) ] for the same parameters. One sees that the region of 
applicability of the leading non-analytic term at this order is considerably smaller, extends 
up to m 2 ~ 0.05 GeV 2 . The relativistic result gives a better description out to larger pion- 
mass values due to a resummation of higher order effects (m* lnm^, m%, etc., terms), which 
ensures the correct analyticity properties. The pion-nucleon sigma-term can be obtained in 
this calculation as cr/v = m^dM^/dm^, taken at physical m w : 

a N = 67-U = 50 [MeV] , (21) 

where the first number refers to the contribution of the low-energy constant Cxjv, while the 
second is the chiral loop correction. 

The dotted curve in Fig. |21 shows the relativistic result for both ttN and it A loops ac- 
cording to Eq. (fTTj) . with slightly re-adjusted low energy constants = 0.870 GeV, 
citv = — 1.1 GeV" 1 . One sees that the 7rA loop gives a substantial contribution for larger 
pion masses and spoils the agreement of our p 3 relativistic calculation with lattice data, 
above m 2 ~ 0.15 GeV 2 . The corresponding sigma-term is 

a N = 85 - 17 - 11 = 57 [MeV] , (22) 

where the numbers refer to the contributions of cjjv, the ttN, and the 7rA loops, respectively. 

In absence of a complete fourth order calculation for both irN and 7rA loop contributions, 
we estimate the higher order terms here in a simple way by allowing for one additional term, 
proportional to m\ in the baryon mass formula as : 

M B = M B 0) - 4 cxb ml + c 2B m\ + chiral loops, (23) 

where the chiral loop contribution is calculated as discussed above in Eq. ()17j) . Such a 
procedure was also proposed before in Ref. |9|, when applying a heavy-baryon formula for 
the non-analytic contribution in the quark mass to lattice results. We see from Fig. |21 
that with this 3 parameter formula, one can obtain a description of the lattice calculation 
with the relativistic tcN + 7rA result up to m 2 ~ 0.5 GeV 2 , using as parameter values : 
= 0.87 GeV, c 1N = -1.1 GeV" 1 , and c 2N = 3.0 GeV" 3 . Although the m\ term only 
contributes about 1 MeV to the nucleon mass for physical pion mass values, its contribution 
at m 2 =0.5 GeV 2 amounts to about 750 MeV. We notice that in Ref. [§| the value of the 
coefficient c 2 n was also found to be large, which signals the importance of further higher- 
order terms. 

Fig. El displays the results for the A mass. As in the nucleon case, the relativistic ttN 
loop result, shown by the solid curve, is able to provide a surprisingly good description of 
the lattice results up to about m 2 ~ 0.4 GeV 2 , using Nlf ] = 1.20 GeV, c 1A = -0.40 GeV" 1 . 
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FIG. 3: (Color online) Pion-mass dependence of the A mass. The solid (black) curve is the 
relativistic ttN loop result, with = 1.20 GeV and c iA = -0.40 GeV -1 . The dashed (blue) 

curve is the leading non-analytic m 3 result (arising from 7rA loops), with = 1.185 GeV and 
cia = —0.75 GeV" 1 . The dotted (green) curve shows the relativistic 7riV + 7rA result for the same 
parameters, whereas upon adding to this result the m% term as in Eq. (j2HJ), with C2A = 2.0 GeV" 3 , 
one obtains the dash-dotted (red) curve. The (red) squares are lattice results from the MILC 
Collaboration 2]. The star represents the physical mass value, which is used in the fits. 

When including the 7rA loops one notices that although the convergence of the relativistic 
calculation (dotted curve in Fig. |3J) is improved in comparison with the leading-nonanalytic 
result (dashed curve in Fig. |3J), its agreement with the lattice results is limited to m 2 ~ 
0.1 GeV 2 . 

As for the nucleon, we estimate remaining fourth order contributions by the form of 
Eq. ()23|) . Using such a three parameter form, the relativistic ttN + 7rA loop calculation 
is able to describe the pion mass dependence of the A mass up to m 2 ~ 0.5 GeV 2 , with 
Mf = 1.185 GeV, c 1A = -0.75 GeV" 1 , and c 2A = 2 GeV" 3 . We note that the coefficient 
C2A is of the same size as for the nucleon and represents a 500 MeV mass contribution to 
the A mass at m 2 = 0.5 GeV 2 . 

In conclusion, we have studied the pion-mass dependence of nucleon and A-isobar masses 
within the framework of a manifestly covariant chiral effective-field theory. We have com- 
puted the one-loop nN and it A graphs and obtained a generic expression for those contri- 
butions to the masses and field renormalization constants, Eq. (J17)) . We were able to obtain 
these generic expressions because of a specific choice of the chiral Lagrangian, where the 
ttNA and 7rAA couplings are constructed to be consistent with the spin degrees of freedom 
counting of the relativistic spin-3/2 field. For the ttNN coupling we adopt a form which is 
similar to the above-mentioned A couplings. The resulting relativistic loop corrections obey 
the chiral power-counting, after renormalizations of the available counter-terms are done. 
The relativistic expressions also contain the nominally higher-order terms, which are nec- 
essary to satisfy the analyticity constraint. As has been shown here on the example of the 
nucleon and A masses, the convergence of the chiral expansion can be improved in this way, 
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without introducing additional free parameters. In particular, we find that the relativistic 
calculation, including only the ttN loops, is able to describe the lattice results for both nu- 
cleon and A masses up to m 2 ~ 0.5 GeV 2 with only one free parameter. Including the 7rA 
loops, however, spoils the agreement with the lattice result. We then estimated the effect of 
the higher-order terms by adding a m\ term to the baryon mass. Using the additional free 
parameter, one is able to obtain a description of the lattice calculation with the relativistic 
tiN + 7i A result up to m 2 ~ 0.5 GeV 2 . While a full fourth order calculation of both ttN and 
7rA loops is a worthwhile topic for future work, the present relativistic chiral-loop calculation 
can be used in the interpolation between full lattice QCD simulations and the experimental 
results. 



APPENDIX A: CONVENTIONS, RULES AND IDENTITIES 

Here we summarize the conventions, Feynman rules, and list a few useful identities used 
throughout this work. 



^ = diag(l, -1,-1,-1), e 0123 = 1, 75 = *7o7i7273, (7 M ) f = tVt , 
ermore, 7's denote Dirac's 7-matrices and their totally-antisymmetric 

- ^[7 ;U ,7 iy ], 7^ a = ^{7^,7°}, 7^ Q/3 = 1 \pi^ va ) y 3 ] . 



Conventions: g , . . , ; , . 

75 = 75. Furthermore, 7's denote Dirac's 7-matrices and their to 
products: 7^ = ^ ua = lil^iT}, 1^ = \[Y ua , 

Propagators: 

SM = (p 2 -ml + ie)-\ 

Sn(p) = (^-M N + te)- 1 = + M N ) (p 2 - M 2 N + ie)" 

Sfip) = -{j> + M A )(p 2 - Ml + iey l vW 2 ^{p). 



• Vertices: 



r. 



(l)o 



■kNN 



(P',P) 



r. 



(l)a,a/3 ( , 



ttAA 



(p',p) 



2M N f n 

h A 
2M A U 

2M a U 



te a ^%p elal5 T a , 



ie al3ea p,p' a T a 



(Ala) 
(Alb) 
(Ale) 



(A2a) 
(A2b) 
(A2c) 



• Identities: 



(A3a) 



ie^ 75 = rf*> = KT^T V - 7 e 7 !/ 7 M 7 <T + 7V7V ~ Yl a Yl e ) , (A3b) 

£ A<ra P X l" £ ^PP P U 7 P = -p 2 (ga/3 + lalfl) + flfiPa + laPpf, (A3c) 

e^ aP x r (g^ - IrY) e MP »Y = -P 2 V% 2 ( P ), (A3d) 

e a ™' a 'p e ,k a , e Pvea p g k a g^ = -g a ^[p 2 k 2 - (p • A;) 2 ] + p a p k 2 + A; a A;V 

-p-k(p a k< 3 +p /3 k a ), (A3e) 

e^p^^e^p^T^ip) = -p 2 k a k?V% 2 (p). (A3f) 
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APPENDIX B: RELATIVISTIC LOOP FUNCTIONS 



Here we explicitly define functions V and W which enter the mass and FRC correction 
formula fT7|) . 

V(jm, 6) = V(jm, 6) - V(0, 6) - V'(0, 6)fi 2 , (Bla) 

V(fi,8) = / dx(R + x) {/j 2 x + (1 -x){R 2 -x)} ln{/i 2 x + (1 - x)(R 2 - x)} 
Jo 

= \{R + a)[(3{^ 2 -2\ 2 )\n^ 2 + a {R 2 -2\ 2 )\n R 2 ~l(a 3 + (3 3 ) (Bib) 

+ 4A 2 + 4A 4 fi(A)] +i/i 4 (ln/i 2 -i)-ii? 4 (lni? 2 -i), 

V(0,6) = ^ + ^R-^R 2 -lR 3 -lR* + lR 5 + ±R 6 (Blc) 

- ±R 5 {R 3 + 2R 2 -2R-6) In R + ±{R 2 - 1) 3 (1 + R) 2 ln|i2 2 -l|, 

V'(0,5) = (d/d/j, 2 )V(fi,5)\^ =0 = -^(7 + 9R + 3R 2 + 9R 3 + 6R 4 ) 

+ i J R 5 (3 + 2i?) InR+^l + lR-'lR 5 -R 6 ) ln|i? 2 -l|, (Bid) 

where R = 1 + 5, f3 = -<5-i(<5 2 -/i 2 ) = ^(l-i^ + zi 2 ), a = l-[3, A 2 = \[5 2 -^ 2 } [(2+<5) 2 -/i 2 ], 
and the elementary function Q is defined as: 



Q(\) = { a 7^' CT \ " U (B2) 
-7=F arctan ^?' A < °- 

Similarly, 

F(/i, 5) = W(ji, 6) - W(0, 5) - W'(0, 8)jj 2 , (B3a) 



= [ dxx {fi 2 x + (1 - x)(R 2 - x)} \n{fi 2 x + (1 - x)(R 2 - x)} 
Jo 

= fa [(3 Gu 2 - 2A 2 ) ln/i 2 + a {R 2 - 2A 2 ) lni? 2 - |(a 3 + f3 3 ) (B3b) 
+ 4A 2 + 4A 4 fi(A)] +i/i 4 (ln y u 2 -l)-ii? 4 (lni? 2 -i), 



- \R & {R 2 - 2R) \nR + ±(R 2 - 1) 3 (1 + R 2 ) In \R 2 - 1| , (B3c) 
W'(0,5) = (d/dti 2 )W(ti,6)\ n= =-j-J7 + 3R 2 + 6R*) 



l/i=0 18 1 

+ f i? 6 \nR + 1(1 - R 6 ) In |i? 2 - 1| , (B3d) 
It is useful to know the expansion of these functions for small /x: 
— zx 4 

V{n,5) = !j[-\nfi + lR(l + 2R 2 )-R 5 \nR+±(l + R 5 ) In \R 2 - 1|] + 0(/i 5 )(B4a) 
W M = g h ^ + W + g) - + i(l ± *) - I* - LP + o0 ,. )(B4b) 

For 5 = 0, the expansion takes a different form: 

47T 5 

^(/x, 0) = —/x 3 + -/x 4 (-i + ln/x) + 0(/x 5 ), (B5a) 

27T 3 

WfaO) = T /x 3 + -/x 4 (-^+ln/x)+0(/x 5 ). (B5b) 
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